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Abstract 

The present work generalizes a nonlocal version of the Polyakov-loop-extended Nambu 
and Jona-Lasinio (PNJL) model to the case of three active quark flavors, with inclusion 
of the axial U(l) anomaly. Gluon dynamics is incorporated through a gluonic background 
field, expressed in terms of the Polyakov loop. The thermodynamics of the nonlocal PNJL 
model accounts for both chiral and deconfinement transitions. Our results obtained in mean- 
field approximation are compared to lattice QCD results for Nf = 2 + 1 quark flavors. 
Additional pionic and kaonic contributions to the pressure are calculated in random phase 
approximation. Finally, this nonlocal 3-flavor PNJL model is applied to the finite density 
region of the QCD phase diagram. It is confirmed that the existence and location of a critical 
point in this phase diagram depend sensitively on the strength of the axial U(l) breaking 
interaction. 

1 Introduction 

Investigating the phase diagram of strongly interacting matter is a persistently challenging theme 
of nuclear and high-energy physics. Ultrarelativistic heavy ion collisions at the Relativistic Heavy 
Ion Collider give evidence that a strongly correlated quark-gluon phase is formed at temperatures 
above 200 MeV, in accordance with lattice QCD calculations (at zero baryon chemical potential) 
which suggest a transition temperature in the same range for chiral symmetry restoration and 
confmement-deconfinement transitions. The ALICE experiment at the Large Hadron Collider 
will shed further light on the region of higher temperatures and low quark chemical potentials. 

Lattice QCD computations can so far not be extended systematically into the region of 
larger quark chemical potentials /x, but expansions around \i = do suggest the existence of a 
critical point at which the chiral crossover transition at small [i turns into a first-order phase 
transition. The precise location of this critical point and possibly even its mere existence are still 
controversial. Given that, at the present exploratory stage, neither lattice QCD nor experiment 
can yet map out the QCD phase diagram over large areas in the (T, jj) plane, it is useful to work 
with models in order to identify the dynamically relevant degrees of freedom. 

Models of the Nambu and Jona-Lasinio (NJL) type [1] have been quite useful for orientation 
in this context [2-12,27-29,41,48] as they properly incorporate the chiral symmetry breaking 
scenario of low-energy QCD. A basic element of such models is the gap equation connecting the 
chiral condensate and the dynamical quark mass, providing a mechanism for spontaneous chiral 
symmetry breaking and the generation of quark quasiparticle masses. Thermodynamic aspects 
of confinement, while absent in the original NJL model, can be implemented by a synthesis with 
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Polyakov-loop dynamics. The resulting Polyakov-Nambu-Jona-Lasinio PNJL model [7-12, 27- 
29, 48] has been remarkably successful in describing the two-flavor thermodynamics of QCD. 
However, this earlier version of the PNJL approach still worked with an artificial momentum 
space cutoff, tInjl ~ (0.6-0.7) GeV, which prohibits establishing connections with well-known 
properties of QCD at higher momentum scales such as the running coupling and momentum- 
dependent quark mass function. Furthermore, thermodynamically consistent results at high tem- 
peratures and densities cannot be achieved using the original (local) PNJL model. In particular, 
a meaningful extrapolation to the high-density region with its variety of color-superconducting 
phases cannot be performed once the quark Fermi momentum becomes comparable to the NJL 
cutoff. The nonlocal PNJL model does not have such a priori limitations. 

In fact the nonlocal two-flavor PNJL model [2] solves this problem by introducing momentum- 
dependent quark interactions that permit realizing the high-momentum interface with QCD and 
Dyson-Schwinger calculations at the level of the quark quasiparticle propagators. The present 
work takes a next major step by extending this nonlocal PNJL model to Nf = 3 flavors, now 
incorporating the strange quark. This step involves a detailed study of the axial U(l) anomaly, 
its role in separating the flavor singlet component of the pseudoscalar meson nonet from the 
Nambu-Goldstone boson sector, and its thermodynamical implications. It will turn out, that 
the nonlocal PNJL model does not suffer from the pathologies mentioned above, and hence it 
does not have any a priori limitations. Therefore, it is well-suited to investigate the high-density 
and high-temperature region of strongly interacting matter. 

From hadron spectroscopy it is well-known that only eight of the nine lightest pseudoscalar 
mesons (the pions, the kaons, and the eta meson) have pseudo-Goldstone boson character. The 
eta-prime meson, on the other hand, has a mass of ~ 958 MeV, about 400 MeV higher than 
the eta mass: the rj' meson is not part of the Nambu-Goldstone multiplet. The reason is the 
axial U(1)a anomaly [30] and its realization in QCD. The flavor singlet axial current, j£{x), is 
not conserved. In fact one has d^jt? oc E a • B a ; i.e. the divergence of the singlet axial current 
is proportional to the product of color electric and magnetic gluon fields. In the interpretation 
given by 't Hooft [31], the axial anomaly is induced by instanton effects and translates into 
an axial U(1)a breaking effective interaction between quarks that has the form of an N{ x N{ 
determinant [32] of right- and left-handed quark bilinears, ± 75)^. For Nf = 3 this is a 
genuine six-point vertex involving all three u, d and s quarks simultaneously. In its local form, 
this 't Hooft determinant interaction has been widely used in previous three-flavor NJL model 
calculations [13-15, 41,48]. Our present work generalizes this U(1)a breaking interaction for its 
use as part of the nonlocal Nf = 3 PNJL model. 

This paper is organized as follows. Sections 2 and 3 introduce and develop the nonlocal 
PNJL model. Its connection with the instanton model is outlined. A comparison with results 
from Dyson-Schwinger calculations of the Landau gauge QCD is performed. Properties of the 
pseudoscalar mesons, such as masses and decay constants, are derived within a systematic 
expansion of the action around the mean-field limit. Fundamental low-energy theorems such 
as the Gell-Mann-Oakes-Renner relation are shown to hold within the nonlocal framework. 
Section 3 proceeds with the thermodynamics of the nonlocal PNJL model using the Matsubara 
formalism. The Polyakov loop <I> is introduced as an (approximate) order parameter for the 
confinement-deconfinement transition. The quarks are coupled to 3> in the usual minimally 
gauge invariant way. The gap equations following from the PNJL effective action determine 
the temperature dependence of the chiral up-, down-, and strange-quark condensates and of 
the Polyakov loop. We also take into account contributions to the pressure beyond mean- 
field approximation, calculating pionic and kaonic quark-antiquark correlations in random phase 
approximation. The nonlocal PNJL model is then applied to the finite density case, introducing 
nonvanishing (quark) chemical potentials. We sketch a simplified version of the QCD phase 
diagram (without inclusion of diquark condensates). The variation of the "critical point" with 
changing strength of the 't Hooft interaction is investigated. Section 4 presents conclusions and 
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an outlook. 



2 Nonlocal three-flavor Nambu— Jona-Lasinio model 

This section extends the previously developed nonlocal NJL model [2] to 7V f = 3 flavors. It is de- 
signed so as to correctly implement the spontaneously broken chiral SU(3)r x SU(3)l symmetry 
of QCD at zero temperature together with the anomalously broken axial U(1)a symmetry. 

2.1 Chirally invariant nonlocal action 

As in our previous work [2] , the construction of the interaction part of the effective quark action 
in Euclidean space 1 is guided by the nonlocal coupling of color currents. A Fierz transformation 
leads to the following generic form of the nonlocal 4-point interaction, 

4 4 t } = £c a /d 4 */ d ' z (* + §) ^ (* - I) G(z) *(*-!) r>(* + f). (2.i) 

The nonlocality distribution, Q(z), is proportional to the gluonic field correlator. As in Ref. [2] 
we restrict ourselves to the diagonal part of the Lorentz tensor representing this correlator. 
The coefficients c a result from a Fierz transform of the color current-current interaction on 
which the action (2.1) is based. The r a are a set of Dirac, flavor and color matrices. Their 
combination is determined by the Fierz transform just mentioned. The quark field is given as 
ip(x) = (u(x),d(x), s(x)) T in terms of the u-, d-, and s-quark fields. Their current quark masses 
are collected in the mass matrix rh q = diag(m u , m^, m s ). 

The coupling term (2.1) of the action shares with full QCD a global chiral U(3)r x U(3)l sym- 
metry. The axial U(1)a subgroup of U(3) x U(3) is broken by the anomaly which will be treated 
by an additional term of the action, to be described in subsection 2.2. Apart from a strictly 
conserved U(l) symmetry associated with baryon number, the remaining chiral SU(3)r x SU(3)l 
symmetry undergoes dynamical (spontaneous) breaking through S-£), down to flavor SU(3) (the 
well-known "eightfold way" ) . Explicit chiral and flavor symmetry breaking corrections are then 
introduced by the quark mass term of the Lagrangian, tprh q ip. 

In Eq. (2.1) we restrict ourselves to the color-singlet pair of scalar and pseudoscalar operators 
as a minimal U(3) x U(3) symmetric combination: 

r ae{o,i,...,8} = (l,i75) x A a , 

with the Gell-Mann matrices {Ao, Ai, . . . , As} specifying a basis in flavor space, where we have 
introduced the singlet matrix Ao := \J\ diag(l, 1, 1). Other less relevant operators (such as 
Lorentz vectors and axial vectors) will be ignored in the present work. (See however the inves- 
tigation of the role of vector current interactions in Ref. [16].) 

Furthermore, we replace the dimensionful nonlocality distribution function Q in Eq. (2.1) by 
a coupling strength G/2 (with mass dimension —2) and a normalized distribution C, i.e. 

Q(z) = |c(z), (2.2) 

with J d A zC(z) = 1. Note that the standard (local) NJL model follows for the limiting case 
C(z) = 8 <yA \z), with the understanding that a momentum space cutoff is introduced in this local 
limit in order to regularize loop integrals. 

x If not stated otherwise, all quantities are expressed in Euclidean space-time, i.e. x = (x4,x) — (ix ,x), 
74 := i7o, etc. 
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2.1.1 Schwinger-Dyson equation 



Following Ref . [2] , the next obvious step is a bosonization of the action 5 int . Such a procedure 
was performed, e.g., in Refs. [17-19] leading to the following Schwinger-Dyson (SD) equation 
(gap equation) in momentum space for the dynamical quark masses M u (p) = Md(p) = M s (p) = 
M(p) in the chiral limit (i. e., m u = = m s = 0): 

"M-^/l^db-rijr^M^. (2.3) 

Here the Euclidean Fourier transform C(p) = J d 4 ze~ ipz C(z) of the distribution C(z) has been 
introduced, with C(p = 0) = 1. The mass function M{p) is interpreted as the momentum- 
dependent dynamical (constituent) quark mass which, in the case of M(p) ^ 0, expresses the 
spontaneous chiral symmetry breaking. The gap equation (2.3) can be solved iteratively relying 
on Banach's fix-point theorem. It turns out indeed that such a procedure leads to a nontrivial 
solution with nonvanishing M (see Fig. 1). 



2.1.2 Four-fermion separable interaction 

The full solution of the SD equation (2.3) is, however, not practical for our present purposes. If 
one were to solve such an integral equation for finite temperature one would be confronted with 
the full complexity of finite temperature SD calculations (cf. Refs. [21-23]). We approximate 
the distribution C by a separable form, replacing 

9M - few ^f/dVc + 

(compare also Ref. [24]). The model using this separable ansatz will be worked out starting from 
the next subsection. A comparison between the full SD calculation and the separable approach 
in Sec. 2.5.1 demonstrates that the separable action is a very good approximation to the full 
SD-formalism for all practical purposes. We can now write the chirally invariant four-fermion 
interaction in a more tractable form: 

5« = -§ J d 4 x [£ (x)j s a (x) + f a {x)j P a {x)] , (2.4) 
where jf and j£ are scalar and pseudoscalar densities, given by 

j s a {x)= U A zi,(x+ Z -) C(z)\^(x- Z -), 

r 1 z (2 - 5) 

£{x) = J d 4 *^(z + |) C(z) i 75 A a V 
The particular functional form of C(z) used in this work will be given in Sec. 2.5 (see also Fig. 2). 



2.2 Axial U(l) anomaly and instantons 

So far, the generalization from the two-flavor to the three- flavor case has been straightforward. 
The more difficult task concerning the extension of the nonlocal NJL approach to three flavors is 
the implementation of the mechanism that breaks the U(1)a symmetry leading to the relatively 
large mass of the rf meson. 

At the level of an effective interaction between quarks, this can be accomplished by introduc- 
ing an interaction term proportional to the Kobayashi-Maskawa-'t Hooft determinant [31,32], 

J d 4 x [det J+(x) +detj-(x)] 
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with 

(J±(x))ij = J dS^(x + |)^(lT75)/C(z)^(x-|), (2.6) 

where JC(z) represents the distribution of the U(1)a breaking interaction strength. For Nf = 3, 
this flavor determinant generates a genuine 3-body interaction (or 6-quark vertex) in which the 
u, d, and s quarks participate simultaneously. 

In his original paper [31], 't Hooft derived the U(1)a breaking interaction starting from 
instantons. The instanton liquid model gives a simple expression for K, as the density of zero 
modes. Its Fourier transform IC(p) = J d 4 ze~ ipz IC(z) is written [34,35] in terms of Bessel 
functions and a characteristic instanton size, d ~ 0.35 fm, as follows: 

1C(P) = *p 2 d 2 ^[l (OKo(0 - h(0Ki(0] with e = ^ . (2.7) 

At this point we can anticipate a result of our present studies (see Fig. 2), namely that the 
distribution C in the nonlocal chiral four-fermion interaction (2.4) turns out to be very close to 
the distribution K. that characterizes the instanton induced interaction (2.6). We can therefore 
set C = K, and work with a universal distribution involving a single scale, e.g., the instanton 
size d. It is then not surprising that the local version of the (classic) NJL model operates with 
a momentum cutoff scale that reflects the inverse instanton size, /Injl ~ 1/d. 

With these preparations and following the steps outlined in Appendix A, the 't Hooft deter- 
minant can be written in terms of the nonlocal scalar and pseudoscalar densities (2.5), leading 
to the following six-fermion part of the action: 

*2 = "f / d ' xA ^ \j*{x)j${x)jS{x) " 3jf(s)#(s)#(s)] , (2.8) 
where the constants A a ^ are expressed in terms of the Gell-Mann matrices according to 

A*/3 7 := ^e ijk e mne {K) irn (\p) jn (A 7 ) w for a, 0, 7 € {0, ... , 8}. 

Here H is the coupling strength of mass dimension —5, representing the six-fermion vertex 
interaction. 



2.3 Three-flavor nonlocal NJL model 

Given the four- and six-fermion couplings in the nonlocal framework, we can now write down 
the Euclidean three- flavor nonlocal NJL action, <Se, that will be the basis of all calculations 
performed in this work. We have 

S E = [ d A x U(x) [-vfdp + rh q ] i;(x) - I [j s a (x)£ (x) + j^(x)j^(x)] + 

H , ( 2 - 9 ) 

-jA*/3 7 [ja {x)j$(x)jS( X ) - 3jZ(x)j£{x)jF(x)] 

where the first term is the kinetic term and rh q = diag(m ?t , m d ,m s ) is the mass matrix with the 
current quark masses m u ,m^,m s ; G and H are constants to be determined and the densities 
JaiJa are given in Eq. (2.5) with the nonlocality distribution C(z) yet to be specified. 

In the remainder of this Section 2 we demonstrate how this approach works in reproducing 
zero-temperature QCD, the nonperturbative vacuum and its lowest quark-antiquark excitations: 
the pseudoscalar meson nonet including decay constants and 77-77' mixing. Thermodynamics and 
the implementation of the Polyakov loop (the step from the nonlocal NJL to the PNJL model) 
will be described in Section 3 and performed as in Ref. [2] by the gauge covariant replacement 
dfi — > d/j + iA^Sfpi, and by adding the Polyakov-loop effective potential to the action. 
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As usual, we start from the partition function 

Z = J &ij^e' SE (2.10) 

and seek to replace the fermionic fields appearing in Eq. (2.10) by bosonic fields. For this reason, 
we introduce 18 bosonic fields a a and 7r a (a £ {0, . . . , 8}) and, additionally, 18 auxiliary fields 
S a ,P a in order to deal with the six-fermion interactions induced by the 'tHooft term. 
Inserting a "one" in terms of delta functions, 

1 = J 9S a ®Pp 5(S a - j s a ) 5{Pp - %) = j ®S a 2>Pp ®a a ^ e Jd 4 za a (s a -^) Jd^^-^ 
the partition function can be written as 

Z = J ®$ 9$ 3>S a @Pp &a a 9np exp j - J d 4 x $ {-vfd» + m q ) $ + a a j s a + irpjf] } 

G H 

(faJa + 3a 3a ) + ~T A ^ O^ltf ~ Zjajpjj) + <T a S a + ITpPf) 



x exp < / d 4 x 



(2.11) 



The second term has expressions quadratic and cubic in the densities j. Applying the inserted 
delta functions dictates a replacement of and j£ by S a and P a , respectively. The first 
exponential, on the other hand, contains terms of the form ipg/tp; hence the path integration 
over the fermionic fields tp and tp can be carried out by standard means, leading to 

Z = J £>a a £>ir a det J J &S a &P a exp | j d 4 x (a a S a + ir a P a ) 

x exp y d 4 x | (S a S a + P a P a ) + ^Aafr (S a S(,S y - 3S a PpPJ I 

where det is the fermion determinant. In momentum space one finds after a simple Fourier 
transformation 

*W) == {pW\p) = H + ™i) ^t^P -p') + C {^Y^j K [°a{p - p) + i75vr Q (p - p')] . 

(2-12) 

Here and in the following we conveniently write C(p) = C(p) for the Fourier transform of the 
distribution C(z). 

2.3.1 Stationary phase approximation 

Owing to the cubic terms in S a and P a the path integration over these fields cannot be car- 
ried out explicitly. Therefore, the stationary phase approximation (SPA) is used, choosing the 
fields S a ,P a so as to minimize the integrand in the bosonized partition function Eq. (2.11). 
Consequently, a necessary condition imposed on the fields is 

3H 

a a + GS a + — Aafo [SpSy - PpP-f] = , 

m (2 - 13) 

TT a + GP a —Aafi-ySpPy = , 

where S a ,P a are now to be considered as (implicit) functions of a a ,7r a . The bosonized action 
can thus be written as 



S£ os = - In det J- J d 4 x [a a S a + 7r a P a + ^ [S a S a + P a P a ] + 

+ -£-Aaf3~ t [S a SpS 1 — 3S a PpPj]^ . 



(2.14) 
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2.3.2 Mean-field approximation, gap equations and chiral condensates 

As a next step we expand the bosonized action, S^ os , in a power series around the expectation 
values of the fields a a ,n a , 



(2.15) 



<r a (x) = a a + Sa a (x), 
ir a (x) = Sir a (x). 

A first constraint is imposed on the scalar fields by charge conservation; i.e., the charge matrix 
Q = diag (|, — |, — |) commutes with the SU(3) generators: [Q, X a ] = 0. This is only possible 
for Ao, A3, As which means in turn that only <to, 03 and as have to be considered (in the isospin 
limit which will be investigated later one has the additional constraint that 03 also vanishes). 
Given these conditions it is useful to introduce 



= diag(a M , a d , a s ) := O A O + 03A3 + 0sAs, 



(2.16) 



and, analogously, S = diag(S u , S d , S s ) = SbAo + S3 A3 + S^Xg. 

Since we have (Tr a ) = (P a ) = to leading order, the action in mean-field approximation 
reads 



E 



Tr In 



p 2 l 3x3 + M 2 (p) 




&iSi + —SiSi j + —S u SdS s 



(2.17) 

where M(p) = diag (M u (p), M d (p), M 3 (p)) with 

M i (p) = m i + a i C(p), (2.18) 

and 13x3 denotes the unity matrix in flavor space and is the four-dimensional Euclidean 
volume. 

The mean-field equations (gap equations) are deduced applying the principle of least action, 

for (7j = <7j (i G {u,d,s}). Taking into account that Si and Pi are both implicit 
functions of <7j, determined through the SPA equations in mean-field approximation (compare 
Eq. (2.13)) one obtains the following set of coupled gap equations 

H 



-GSi 
-8N C 



~£-ijk£-ijk Sj Sfc 



Mi(p) 



p 2 + Mf{p)' 



(2.19a) 
(2.19b) 



Finally, the chiral condensate (qq) can be calculated from 5g os using the Feynman-Hellmann 
theorem, by differentiation with respect to the current quark mass m q . Equivalently, one can 
use the definition 



(qq) 



-iTr lim 



S F (x,y)-SP(x,y) 



with the full fermion Green function, 



S F (x,y) = J 



d 4 p 
(2vr) 4 



jp{x-y) 



Sf(p) , 



and the free quark propagator S F °^ subtracted. Using 

1 



S F 



this leads to 



(qq) = -4N C j 



-f + M q (p) 

d 4 P r M q (p) 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



p 2 + Mq (p) p 2 + m 2 

Note that M(p) — > m q for large p. The subtraction makes sure that no perturbative artifacts 
are left over in (qq) for m q / 0. 
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2.3.3 Second-order corrections and meson masses 

In order to calculate the masses of the pseudoscalar mesons, we go beyond mean-field approxima- 
tion and consider second-order corrections to the mean-field action, extracted from a functional 
Taylor expansion, 

5 E } = \ J d 4 xd 4 i/ ^ 5a a (x) 5ap(y) + ^ J d 4 xd 4 ?/ ^ 5-K a {x) 5irp(y), 

where the second derivatives, Scr ^ etc., are evaluated at the mean- field values a a (x) = a a , 

etc. We now focus on pseudoscalar mesonic excitations and change the basis according to 
TTij = (^ a ^a)ij- This gives a standard representation of the pseudoscalar meson octet: 



(*)« 



K~ 



7T 



_2rj8 I _W 

V6 ^ V~3 



(2.24) 



Defining analogously a matrix a for the scalar mesons, the fermion determinant, Eq. (2.12), can 
be written as 



J(p, P ') = {-f + m q ) (2ir) 4 5^(p - P ')+C 



p — p 



V2[a(p- P ')+i l5 7t(p- P ')]. (2.12') 



Next, we calculate the derivatives appearing in the Taylor expansion. Some caveats of the 
calculation are outlined in Appendix B. The resulting second-order contributions to the action 
are given by 



:(2) 



with 



where 



d 4 p 
(2vr) 4 



G t jM {p) 5(T io(P) S °ke(-P) + G ijM (p) STTijip) 5ir M (-p) 



G f jM (p) = U % 6ii 5 i k + ( r tj,kt) > 



-8N C 



d 4 g 
(2tt)' 



c 2 ( q ) 



g+ • q~ T MM + )M 3 (q-) 

[q+ 2 + M 2 {q+)] [q -2 + M 2 {g - ) y 



q ± = q ± | , and (r^) 1 is defined as the solution of the system 



(2.25) 

(2.26) 
(2.27) 



ij,kt 



(2.28) 



The meson masses can now be determined by writing the second-order term of the action, 
Eq. (2.25), in the physical basis as 



:(2) 



1 



d 4 p 



{gap 2 ) [*°(p) A-p) + ^ + (p) (-p)] + 



/ 2 J (2tt) 4 

+ G K {p 2 ) l2K°(p)K\-p) + 2K+{p)K-(-p)] + 

+ G 88 (p 2 ) m(p) m(-p) + G 00 {p 2 ) T]o(p) rjo(-p) + 2G 08 (p 2 ) Vo(p) Vs(-p)} , 
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where the functions Gp are defined according to Eq. (2.26). If one considers only the isospin 
symmetric case, m u = rrid, then one has 



GAP 2 ) 
G K (P 2 ) 
G 88 ( P 2 ) 
G 00 (P 2 ) 
G 08 ( P 2 ) 



G + §S S \ r //- (/ r 



1 

3 
1 

3 
V2 
3 



6G - AHS U - 2HS S 
2G 2 - GHS S - H 2 Sl 

6G + 4HS U - HS S 
2G 2 - GHS 8 - H 2 S 2 
H(S S — S u ) 



2G 2 - GHS* - H 2 S 2 



+ n- u (p 2 ) + zn; s (p 2 ) 
+ 2nz u {p 2 ) + n- s { P 2 ) 
+ n- u ( P 2 ) - n- s { P 2 



From the construction of the action it is clear that the functions Gp correspond to the inverse 
pseudoscalar meson propagators. The corresponding masses are given by the poles of these 
propagators or, equivalently 



G P (-m%) = 0, for P £ir,K, rj. 



(2.29) 



Finally, we perform a last basis change in order to produce the physical r] and rf mesons. We 
introduce the mixing angle 6 = 8(p 2 ) and write 



?? = ?? 8 cos d v - % sin V v , 
rf = rjs sin 9 ri i + rjo cos 9^/ 



(2.30) 



where 6 V = 6(—m 2 ),6 ri / = 9(—m 2 ,). Introducing the (inverse) rj and rj' propagators G v and Gy, 
respectively, instead of GotbGss^os we obtain for the mixing angle 



tan 29 (p 2 ) = 



2G 08 (p 2 



G 00 (P 2 ) - G 88 (p 2 



(2.31) 



and therefore: 



G V (P 2 



G v ,(p 2 



G 88 + Gi 



oo 



G88 + Gi 



00 



+ \jGl 8 +y 



( Gqq — G 88 
2 



(2.32) 
(2.33) 



2.3.4 Renormalization constants 

One usually introduces renormalized fields 2 , (p{p) = Z^ 1 ^ 2 ip(p), so that the quadratic part of 
the Lagrangian can be written as 



1 



jSfW = -(p? + TO }) <p[p)<p{-p), 

and it is directly evident that masses are identified with poles of the propagators. From this it 
is easy to obtain an explicit expression for the renormalization constants, namely 



dG P (p 2 



dp 2 



, for P G 7T, K, rj. 



(2.34) 



2 Here <p(p) stands generically for any of the fields n a (p), . 
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2.3.5 Decay constants 

In this subsection we calculate the (pseudoscalar) meson decay constants defined as 

(0\Ji a (0)\Mp)) = if*HP» ^ <0|^, Q (0)Mp)) = ifa(3Zl /2 Pll , (2.35) 

where a (x) = ^(:r)7 M 75 4f denotes the axial-vector current. In Appendix C of Ref. [2] we 
outline in some detail the calculation of the (unrenormalized) matrix element (0| a (0)\Tr^{p)). 

In order to calculate this matrix element one has to gauge the nonlocal action in Eq. (2.9). 
This requires not only the replacement of the partial derivative by a covariant derivative, 

where A% (a € {0 ... ,8}) are a set of axial gauge fields, but also the connection of nonlocal 
terms through a parallel transport with a Wilson line, 

W(x, y) = V exp j ^ J da 75 \ a A%(x + (y - x)a) (y M - x^) j , 

where we have chosen a straight line that connects the points x and y. This means that expres- 
sions of the form ip(x)0(z)ip(y) (where 0{z) is an arbitrary field operator) in the action <Se, 
Eq. (2.9), have to be replaced by ip{x) W(x, z) 0{z) W(z, y) ip(y), which guarantees the (local) 
gauge invariance of the underlying Lagrangian. It turns out that the only term that is eventu- 
ally affected by the gauging is the fermion determinant si , Eq. (2.12), which then becomes, in 
coordinate space, 



: [~ [ $y + ^5KA a + m q j 8(x - y)+ 

/ \ / \ / \ ( 2 - 36 ) 

+ C(x-y)W x, X -p-) r a J X -±l) W X -±l,y 



(x,y) 



Here r a stands for either r a = X a or r a = i75A Q , and ip a accordingly for either a scalar field, 
a a , or a pseudoscalar field, ir a . 

The desired matrix element then follows from the gauged fermion determinant according to 



1 1 11 t . 1 . . . S 2 In det srf G 

<0|J3U(0)hr*<p» = - SMp)SA « {t) 



(2.37) 

t 

After a lengthy evaluation of the functional derivatives, we obtain 



A=0 

t=o 



+ 



d? 2 qf + M?(q+] 



xti< dag M C[q--a) 



^ V 2 7 (^ 2 +M|(^))(^ 2 +M2 fe )) J ' 



(2.38) 



with 



P P 
q£ = Q+~0- -a), q- = q--(l + a) 

I I ( 2 -39) 
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Now, the decay constants can be derived from the expression (2.38) and their definitions, 
Eq. (2.35), by contraction with p^, hence 

z -i/2 

U = iPM<OK«(0)l^(p)>-fb-> ( 2 ' 4 °) 

evaluated at the corresponding mass p 2 = —mi. Owing to the properties of the Gell-Mann 
matrices one has / Q/ g = Sapf-w for a G {1,2,3} and f a p = SapfK for a G {4,5,6,7}. On the 
other hand, for the 0- and 8 component we obtain 

h,{p 2 ) = \[2fss{p 2 ) + fuu{p 2 )] 
f00(p 2 ) = ^[2fuu(p 2 ) + fss(p 2 )] 

MP 2 ) = f 80 (p 2 ) = [fuu{p 2 ) ~ fss(p 2 )] ■ 

2.4 Chiral low-energy theorems 

In this section we derive the Goldberger-Treiman and Gell-Mann-Oakes-Renner relations ex- 
plicitly from the nonlocal NJL model presented in this work. For this aim, we expand the meson 
self-energy contribution II~ U , Eq. (2.27), up to first order in the current quark mass m u and the 
momentum p 2 : 

n- u (p 2 ,m u ) = ^ - 2 -^m u + Z-y. (2.41) 

The first term on the right-hand side follows immediately from Eq. (2.27) by setting m u = p 2 = 
and using the gap equation (2.19b) in the chiral limit, m u = (index "0"). The second term can 
be recovered by writing C(p) = J- (M u (p) — m u ) in Eq. (2.27) and using the definition of the 
chiral condensate, Eq. (2.23). The last term follows from the definition of the renormalization 
constant Z n , Eq. (2.34). 

Next we notice by expanding the pion decay constant, Eq. (2.38), that only the term in the 
first line of this equation contributes to order 0(p 2 ). One finds 

hm f uu (p 2 ) = -a ufi Z-} ) . (2.42) 

Using Eq. (2.40) this implies 

fw,0 = VufiKfi > ( 2 - 43 ) 

which is nothing but the Goldberger-Treiman relation. Finally, multiplying both sides of the 
expansion (2.41) with G+y and identifying the pion-mass definition, Eq. (2.29), on the left-hand 
side, and the gap equation (2.19a), on the right-hand side, we get 

f 2 ,o m l = -m u (uu + dd} , (2.44) 

which is the Gell-Mann-Oakes-Renner relation. 

We have thus demonstrated that our nonlocal NJL model is consistent, as expected, with 
fundamental low-energy theorems based on chiral symmetry. At the same time the nonlocal 
model reduces to the local NJL model results when C(p) is chosen as a theta function. 
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p [GeV] 



Figure 1: Comparison of the self-consistent 
solution of the Schwinger-Dyson equation 
(2.3) (solid line) to the solution of the 
gap equations (2.19) derived for a separa- 
ble interaction (dashed lines) for values of 
the chiral condensate (uu) ranging between 
-(260 MeV) 3 and -(280MeV) 3 (gray band). 




0.5 1.0 1.5 2.0 2.5 
P [GeV] 



Figure 2: Comparison of the distribution C(p) 
used in the nonlocal NJL model (solid line) 
and the one derived from the instanton model 
with instanton size d = 0.35 fm (dashed line). 
The dash-dotted line shows the function C(p) 
in a local NJL model with Euclidean four- 
momentum cutoff. 



2.5 Parameter fixing and numerical results 

Now that all relevant formulas have been derived in the preceding subsections, we fix the model 
parameters in order to quantitatively reproduce physical observables. Apart from the coupling 
strengths G and H of the four- and six-fermion interaction vertices, this involves a specification 
of the nonlocality distribution C. In preparation of this input we are guided by results of 
Schwinger-Dyson QCD calculations and return briefly to the issue of the separable ansatz for C 
used in the present approach. 

2.5.1 Schwinger-Dyson calculations vs separable approximation 

For simplicity we restrict ourselves in this subsection to the case of two degenerate (up- and 
down-quark) flavors and H = 0. The full SD equation is then given by Eq. (2.3), while the 
separable equation follows by replacing C(q — p) — > C(q) ■ C{p). In this case, the mass equation 
can be written as M(q) = m u + a u C(q) (see Eqs. (2.17) and (2.19a)), with m u = the current 
quark masses and a u = ad the mean-field values of the a fields. Furthermore, for this test case 
and exclusively in this subsection, we choose a Gaussian for the momentum distribution function 
C in both the full SD expression and the separable ansatz. The parameters G, m u and the width 
of the Gaussian are fixed such as to reproduce the pion mass, m n = 140 MeV, the pion decay 
constant 3 , f n = 92 MeV and the chiral condensates ranging between (uu) = (dd) ~ —(260 MeV) 3 
and -(280 MeV) 3 . 

Figure 1 shows results from the full SD and separable nonlocal NJL approaches in com- 
parison. It is evident that the results coincide at the five-percent level. One should recall at 
this point that the volume element in four dimensions is d A p ~ p 3 dp, so that details of the 
low momentum behavior of C(p) in the integrand do not dramatically influence the value of the 
integral in the gap equation. 

3 For the SD equation we use the formulae given in Ref. [17], while for the separable ansatz we use formulae 
(2.29) and (2.40). 
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G 


H l m u 


m s 


(0.96 fm) 2 


-(0.63 fm) 5 3.0 MeV 


70 MeV 



Table 1: Scenario I parameter set of the Nf = 3 nonlocal NJL model. 



(uu) = (dd) (ss) 


M u = M d 


M s 


-(0.282 GeV) 3 


-(0.303 GeV) 3 


362 MeV 


575 MeV 





m K 


m v 


mw 


u 


fx 


6 V 




138 MeV 


487 MeV 


537 MeV 


954 MeV 


83.4 MeV 


104.1 MeV 


3.3° 


-29.1° 



Table 2: Calculated physical quantities using the scenario I parameters of Table 1. 



2.5.2 Parameters of the nonlocal 3-flavor NJL model 

Apart from the nonlocality distribution C(p) with its characteristic scale, the parameters to be 
fixed are the coupling strengths G and H and the current quark masses m u (= ra^) and m s . 
Consider first the distribution C(p) again. Its asymptotic form is constrained by QCD through 
the high- momentum behavior of the quark mass function M(p) to leading order in the operator 
product expansion: 

M(p) oc -^p-(uu + dd) atp>/l~lGeV. 
p z 



-l 

We use the leading order expression a s (p 2 ) = 4ir /3nln-r^ — which is sufficient for our pur- 

L a qcdJ 

poses. We set (3q = 9 (with three active flavors) and Aqcd = 0.25 GeV which reproduces 
a s = 0.12 at mz = 91.2 GeV. At smaller momentum scales we are guided by results from 
Schwinger-Dyson calculations (in Landau gauge) and by extrapolations of lattice QCD data for 
M(p) [37]. 

An alternative choice is to use an instanton model. All these options lead to very similar 
shapes of C{p). For the actual calculations we adopt the parametrization introduced in Ref. [2]: 

for p 2 < A 2 

<V) = t «s(p 2 ) f 2 . A2 (2-45) 

const. • — for p z > A z . 
\ p z 

The normalization condition C(p = 0) = 1 and the matching at p = A fix the remaining constant 
once d is chosen. Setting d = 0.35 fm guided by a typical instanton size, C(p) closely resembles 
the instanton based distribution (2.7), as demonstrated in Fig. 2. The resulting matching scale 
in Eq. (2.45) is A = 0.85 GeV. 

Having determined C (p) we proceed by reproducing masses and decay constants of the pseu- 
doscalar meson nonet. Two scenarios with marginally different coupling strengths G (leaving 
C(p) and all remaining parameters unchanged) will be considered. "Scenario I" optimizes the rj' 
sector including rj-r/' mixing. "Scenario II" provides a best fit to pseudoscalar octet observables. 

Choosing the parameters of scenario I as given in Table 1, one finds the values of the pseu- 
doscalar masses 4 , decay constants and ^'-mixing angle as shown in Table 2. The current quark 
masses are consistent with those listed in the PDG table [38] at a renormalization scale of about 
2 GeV. The rj mass is very close to its experimental value. The same is true for the ratio of 



4 Note, that the uw-threshold is lower than the r\ mass. Hence, the integrals determining the r;'-mass might be 
ill-defined owing to poles in the integration region. Therefore, in fixing the ?7'-mass, we apply the regularization 
method described in Refs. [3,33]. 



13 



G 


H m u 


m s 


(1.04 fm) 2 


-(0.63 fm) 5 3.0 MeV 


70 MeV 



Table 3: Scenario II parameter set of Nf = 3 nonlocal NJL model. 



(uu) = (dd) 


(88) 


M u = M d 


M s 


-(0.304 GeV) 3 -(0.323 GeV) 3 


468 MeV 


694 MeV 





m K 


m v 




U 


!k 


6 V 


6 v r 


139 MeV 


495 MeV 


547 MeV 


964 MeV 


92.8 MeV 


110.1 MeV 1.9° 


-22.3° 



Table 4: Calculated physical quantities using the parameters of Table 3. (These values together 
with the parameters of Table 3 are referred to as "scenario II".) 



the decay constants fx/f-K = 1-25 (compared to the experimental (fx/ fn)exp = 1-22). The pion 
decay constant /„-, though, is approximately 10 % off its experimental value but close to its value 
at the chiral limit. 

Furthermore, it is instructive to compare our result for the ry-ry'-mixing angle, 9 r f = —29.1°, 
to the empirical value. The most recent analysis [36] gives 5 9 = —29.0° and agrees perfectly 
with our result. Note, however, that in Ref. [36] contributions from the gluon condensate are 
included which our model does not explicitly account for. The left part of Fig. 3 shows the 
momentum dependence of the resulting dynamical up-quark mass, M u (p), compared to lattice 
data from Ref. [37]. 

The parameters of scenario II (Table 3) differ from those of scenario I only by a four- 
fermion coupling constant G that is about 15 % larger. We see from the calculated quantities 
given in Table 4 that the pseudoscalar meson masses and decay constants now agree perfectly 
with their empirical values. On the other hand, the magnitudes of the chiral condensates and 
the dynamical quark masses M„(0) and M s (0) increase, making them less compatible with 
common phenomenology. This can easily be understood from the Gell-Mann-Oakes-Renner 
relation (2.44) recalling that the current quark and pion masses for scenario I are the same 
as for scenario II, while the value of the pion decay constant of scenario II is increased. The 
momentum dependence of the dynamical quark mass, M(p), is shown on the right-hand side of 
Fig. 3. In particular, the rj-rj' -m.ix.ing angle 0^ = —22.3° now differs by 20 % from the deduced 
empirical value in Ref. [36]. 

Finally, a comparison with the standard local NJL model is instructive. From Refs. [4, 13-15] 
one finds the gap equations 

H 

M u = m u - G(uu) - — (uu)(ss) 
H 

M s = m s - G(ss) - —(uu) 2 . 

The equivalent coupling strengths G and H of the local model can be evaluated by comparison 
with Eq. (2.19a). One derives for scenario I: G = GtI\ ps HGeV" 2 and H = H v « 

400GeV~ 5 , values that lie well in the ballpark of typical local approaches [4,9-15]. 

We will comment further on the influence of the different model parameters in Sec. 3.4 when 
dealing with thermodynamics resulting from the nonlocal approach. 

5 Note the different definitions of the ?7-?7'-mixing angle in this work and in Ref. [36]. The cited number 
9 = —29.0° has, however, already been translated to the definition, Eq. (2.30), of the mixing angle used in the 
present work. 
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Figure 3: Momentum dependence of the dynamical (constituent) up-quark mass M u compared 
with lattice data extrapolated to the chiral limit (from Ref. [37]). The left figure shows the 
result for scenario I, the right figure the result for scenario II. 

3 Thermodynamics of the nonlocal PNJL model 

In this section we apply the nonlocal NJL scheme in the modeling of QCD thermodynamics. 
Temperature T and quark chemical potential /U are introduced using the Matsubara formal- 
ism [43]. Important features of the confinement-deconfinement transition are introduced by the 
Polyakov loop which serves as an order parameter for confinement in the pure gauge case [39, 40]. 
Coupling the NJL quarks to a gluonic background field, expressed in terms of the Polyakov loop, 
allows describing both the chiral and the deconfinement transitions simultaneously. The frame- 
work for this is the Polyakov-loop-extended nonlocal Nambu-Jona-Lasinio (nonlocal PNJL) 
model. 

3.1 Polyakov loop and its effective potential 

The renormalized Polyakov loop, ($), is defined through 



where tr c denotes the trace over color only and L is the Polyakov loop, a timelike Wilson 
line connecting two points at t = and t = —i(3 in imaginary time, with periodic boundary 
conditions: 



Here V is the path-ordering operator and A4 := A^t a = iAgt a is the fourth component (in 
Euclidean space) of the gluon field 6 , and [5 = 1/T. 

As shown e. g. in Ref. [39], (<fr) = implies an infinite free energy, corresponding to confine- 
ment, while ($>) = 1 implies a vanishing free energy of the system and, thus, deconfinement. Ac- 
cording to Ref. [40], (<£) strictly serves as an order parameter for the confinement-deconfinement 
phase transition only in absence of quarks. But even when quarks are present (<3?) is still useful 
as an indicator for a rapid crossover transition. 

Without loss of generality we can restrict ourselves to the diagonal elements of the SU(3) Lie 
algebra (the so-called Polyakov gauge). We neglect spatial fluctuations and treat A\ and as 
constant fields, so that the Wilson line (3.6) simplifies to a phase factor e 1 ^^ ■e 1 ( X4 ~ y4 ^ A 4 t3+A 4 t8 \ 

6 The color SU(3) matrices t a ,a G {1, . . . , 8} are defined through the Gell-Mann matrices A Q , a € {1, . . . , 8} as 



(Hx)) = —(tT c [L(x)}} 



(3.1) 




(3.2) 
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a 


ai 


a2 


03 


6 4 


3.51 


-2.56 


15.2 


-0.62 


-1.68 



Table 5: Parameters of the Polyakov potential U (from Ref. [10]). 



The contribution of the nondiagonal elements of the SU(3) Lie algebra can be included by 
carrying out the group integration, leading to the following Haar volume 

J(<h,<h) = y— — / J] d & = ^2 ( cos (^3) " cos(^ 8 )) 2 sin 2 (0 3 ) , (3.3) 

SU ( 3 ) 7 ie{l,2,4,5,6,7} 

A 3 ' 8 

where we have set fog = /?—§—. This volume can be written in terms of the Polyakov loop $ 
and its conjugate 



1 - 6$*$ + 4 (V 3 + $ 3 ) - 3 (<T$) 2 



(3.4) 



with $ = ^tr c [exp (i((/>3A3 + (foAs))]- This procedure leads to the construction of the following 
effective potential U [7, 8] that incorporates the effects of the six nondiagonal gluon fields: 



The coefficients are parametrized as 



1 - 6$*$ + 4 (V 3 + $ 3 ) - 3($*$) 2 



(3.5) 



6 4 (T) = 6 4 



The first term on the right-hand side is reminiscent of a Ginzburg-Landau ansatz. The values 
of the coefficients are taken from Ref. [10] and listed in Table 5. The potential U is mani- 
festly invariant under transformations with elements of Z(3), the center of SU(3), which is the 
underlying symmetry behind the confinement-deconfinement phase transition. 

The parametrization (3.5) of the Polyakov-loop effective potential U is applicable at temper- 
atures T up to about twice the critical T c . At higher temperatures, transverse gluon degrees of 
freedom - not covered by the Polyakov loop - begin to be important. Other parametrizations of 
U are possible. An example is the two-parameter ansatz [41] based on the strong coupling limit. 
These two versions of U differ at high temperatures but produce very similar pressure profiles 
[41] at T < 2T C , the temperature region of primary interest in the present study. 



3.2 Coupling of quarks and Polyakov loop 

The coupling of the quarks and the Polyakov loop <E> is introduced by the minimal gauge coupling 
procedure applied to nonlocal field theories. This implies, first, a replacement of the partial 
derivative <9 M by a covariant derivative = <9 M — iA^, or in momentum space — > p M + A^. 8 
Furthermore, in close analogy to the gauging of theories on a discrete lattice (see e. g. Ref. [42]), 
one introduces the gluonic fields writing a Wilson line (compare Sec. 2.3.5) 



Wa(x,v) 



(3.6) 



7 From here onward we omit angled brackets for notational simplicity, i. e. we write (<J>) 
8 The coupling strength g is absorbed in the definition of the fields A^. 
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between the (nonlocal) fermionic bilinears, i.e. ip(x)ifj(y) — > ip(x)yVA(x,y)if)(y). 
In the present context we set 

A„ = 5^{A%+A\t S ), 
3 8 

following our previous discussion, with constant fields A 4 ' . The gauge invariant replacement in 
quark momentum space is then simply 



Pa -»• Pi - {A\U + A 8 4 t 8 ) 



(3.7) 



keeping the three-momentum p unchanged. 

The next step is now the application of the Matsubara formalism as described in the standard 
literature (e.g., Ref. [43]). The variable p^ is replaced by the fermionic Matsubara frequency 
oj n = (2n+l)-7rT, n G Z and taking into account the gauging (3.7). The resulting thermodynamic 
potential in mean-field approximation is 9 



-1/ 



£ (*fSf + ^S f S f ^j + *S u S d S s I + W(<&, <&*, T) , 



(3.8) 



I fe{u,d,s} 



with 

s _1 (iw n ,p; 



iw„7o -7-p-M-i(A 4 + i/i)7o 

iw„7o — 7 ■ p — M* + i(A 4 + i/t)7o 



(3.9) 



where the momentum-dependent dynamical mass matrix M is diagonal in color and flavor space, 



M 



fdmg c (M(u- n ,p), M{cj+ n ,p), M(w° \ 
diag c (M(a;- n ,p),M(a;+ n ,p),M(a;O n ,p)) 

diag c (M(^ n ,#),M(a;+ n ,p),MK n ,p))y 



with = ta 



i/U/± A|/2 — A 4 /(2\/3), = w„ — i/x/ + A|/\/3. Note, that a quark chemical 
potential ft = diag f (/i u , /i^, fi s ) has been introduced that will, however, become important in 
Sect. 3. 6. 10 The trace may be further simplified leading to 



^ = -* T E EE 

fe{u,d,s} i=V>,±nd1, 

1 



d?p 
(2^)3 



Re On 



E 

f£{u,d,s} 



VfSf + ^S f Sf\ + ^S u S d S s + W(d>, d>*,T) . 



(3.8') 



This is the thermodynamic potential of the nonlocal PNJL model in mean-field approximation. 
The auxiliary scalar fields Sf are determined by the SPA conditions, Eq. (2.19a). 



9 Note an extra factor \ because of the doubling of the degrees of freedom in Nambu-Gor'kov space. 
10 The chemical potential enters the 4-component of the argument of the mass function according to the same 
reasoning applied above for the gauging of the model, taking into account that fi can be considered as an imaginary 
potential in Euclidean space. 
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3.3 Gap equations in mean-field approximation 

Once the thermodynamic potential is calculated, the fields a u = ad,cr s and ^4,^4 can be 
determined by requiring thermodynamic potential to be stationary. The necessary conditions 
are given by the gap equations 

dtt _ an _ dQ _ dQ 

d^-dFs~dAZ~dA*~ 0, (3 ' 10) 

together with the stationary phase approximation equations (2.19a). First, we limit ourselves to 
the zero-density, i. e. jl = case. Following Refs. [9-11] we have then = <I>* in the mean-field 
approximation and, consequently, A\ = 0. 

Fig. 4 shows the results for the temperature dependence of the chiral up- and strange-quark 
condensate and of the Polyakov-loop using the parameters given in Table 1 (scenario I) and 
Table 5. This figure illustrates once more, as already demonstrated in Refs. [2,10,12], the 
entanglement of chiral dynamics and Polyakov loop degrees of freedom, a characteristic feature 
of the PNJL approach. In the absence of a coupling between quark quasiparticles and Polyakov 
loop the chiral transition (for Nf = 3 flavors) and the first-order deconfinement transition (of 
pure gauge QCD) appear at very different critical temperatures (T c hi ra i ~ 110 MeV for the chiral 
transition and To ~ 270 MeV for deconfinement). The presence of quarks breaks the Z(3) 
symmetry explicitly and turns the first-order deconfinement phase transition into a continuous 
crossover. The quark coupling to the Polyakov loop moves this transition to lower temperature. 
At the same time the chiral transition (with explicit symmetry breaking by nonzero quark mass) 
turns into a crossover at an upward-shifted temperature, just so that both transitions nearly 
coincide at a common temperature T c ps 200 MeV. 

This symmetry breaking pattern seems also to be realized in recent lattice QCD results using 
staggered fermions [20] where a common chiral and deconfinement transition temperature T c = 
(196 ± 3) MeV is observed. Fig. 4 shows these lattice data for orientation. In the latest work of 
the same collaboration [?], domain- wall fermions are used instead of staggered fermions, leading 
to T c = 171(10)(17) MeV, a value which is consistent with both T c = 196 MeV and alternative 
lattice computations [44] that find a lower chiral transition temperature T c ~ 150 MeV and a 
displacement from the deconfinement transition. 

Compared to the Nf = 2 flavor case (Ref. [2]) we observe only minor changes at this point, in 
particular the transition temperature decreases slightly from T c 2f = 207 MeV to T c 3f = 200 MeV. 
The transition temperature T c can be decreased further if the response from quark effects is 
included in the Polyakov-loop effective potential, leading to a lower To in its parametrization. 
According to Ref. [?] one has T = 190 MeV for 2 + 1 flavors instead of T = 270 MeV for the 
pure gluon case. 

3.4 Parameter dependence 

The chiral and deconfinement transition pattern shown in Fig. 4 changes only marginally when 
scenario I (with parameter set listed in Table 1) is replaced by scenario II with a slightly 
larger coupling G. It is instructive also to examine the dependence on other parameters such 
as the current quark mass m u . In the chiral limit, m u = — > 0, the chiral condensate 
displays a second-order phase transition as expected. Explicit chiral symmetry breaking with 
fnu,d 7^ turns this into a crossover transition as evident from Fig. 4 for m u ^ = 3 MeV. 
Increasing the quark mass to m u ^ = 10 MeV makes the crossover softer at T > T c while 
leaving the condensate unaltered at temperatures below T c . The reason is that, below the 
transition temperature, the dynamical quark mass M(p) entering the chiral condensate in the 
finite temperature generalization of Eq. (2.23) is dominated by the large scalar field a u ~ 
0.4 GeV. Changes of the light quark mass m U) d within a 10 MeV range are not important 
compared to that scale, whereas they become more prominent above T c where a u drops rapidly. 
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Figure 4: Results of nonlocal PNJL calculations of the chiral and deconfinement transition 
pattern using the parameter set of scenario I. Left solid curve: temperature dependence of the 
chiral condensate (uu) = (dd). The strange-quark condensate (ss) is shown as the dashed curve. 
Right solid curve: temperature dependence of the Polyakov loop <E>. Left dotted curve: light 
quark condensate without coupling to Polyakov loop. Right dotted curve: Polyakov loop in the 
absence of quarks (pure gauge QCD). Also shown are lattice results for the chiral condensate 
and the Polyakov loop from the "hotQCD" collaboration [20]. The temperature is given in units 
of the transition temperature T c = 200 MeV. 

The softening of the strange-quark condensate (ss) above T c , as seen in Fig. 4, is much more 
pronounced, given the larger s-quark mass m s ~ 70 MeV. 

Corrections to the behavior of the chiral condensate and the pressure below T c come primarily 
from thermal pions (and kaons) as will be discussed in the following subsection. 

3.5 Beyond mean field: mesonic corrections 

So far the calculations have been performed in the mean-field approximation in which the pres- 
sure P = — Vt is determined by the quarks moving as quasiparticles in the background provided 
by the expectation values of the sigma fields, a u = <Jd and a s , and of the Polyakov loop <F In or- 
der to get a realistic description of the hadronic phase (at temperatures T < T c ), it is important 
to include mesonic quark-antiquark excitations. The hadronic phase in the absence of baryons is 
dominated by (the light) pseudoscalar mesons (pions and kaons). The quark-antiquark contin- 
uum is suppressed by confinement. However, as pointed out in Ref. [46], mesons described within 
the standard (local) PNJL model can still undergo unphysical decays into the quark-antiquark 
continuum even below T c . In the nonlocal PNJL model, such unphysical decays do not appear 
by virtue of the momentum-dependent dynamical quark mass 11 . This means that the pressure 
below T c is basically generated by the pion and kaon poles of the corresponding one-loop qq 
Green functions, with their almost temperature independent position. Therefore the calculated 
pressure below T c corresponds to that of a boson gas with constant masses. 

To include the mesonic contributions to the pressure in our nonlocal PNJL model we can 

11 With the possible exception of the 7/'-meson which will not be considered in this section. 
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T [GeV] 

Figure 5: Pressure (in units of T 4 ) calculated in the nonlocal PNJL model as a function of 
temperature on an absolute temperature scale. Solid curve: full calculation (i. e. mean-field 
result plus mesonic corrections). Dash-dotted curve: mean-field result (no mesonic corrections). 
Dashed curve: mean-field plus pionic and corresponding scalar modes. 



basically use the formalism described in Ref. [46]. One has to calculate Gps,s(i / m,P ), Eq. (2.26), 
and, in particular, the quark loop contribution to the pseudoscalar (PS) and scalar (S) mesonic 
self-energies -/Ips,s( z/ m;P ) (where v m = 2TrmT,m £ Z is the bosonic Matsubara frequency and 
p is the momentum of the incoming meson), given in Eq. (2.27), at finite temperature. Using 

the replacement J — > T^ nGZ J and the rules (3.7) this can be carried out easily. 
The additional contribution of mesonic quark-antiquark modes to the pressure is given by a ring 
sum of random phase approximation (RPA) chains, investigated in Ref. [47] and leading to the 
expression 

P meS on(T) = -T £ [-^HGM(v m ,p)}, (3-11) 

M=PS,S me% ^ ' 

where du is the mesonic degeneracy factor [d,M — 3 for pionic and du = 4 for kaonic modes). 
Because of the momentum dependence of the nonlocality distribution C(p) and the dynami- 
cal quark masses M q (p), integrations and summations in Eq. (3.11) can only be carried out 
numerically. 

Results for the pressure in the presence of pion, kaon and scalar modes are presented in 
Fig. 5. 12 Apart from the full result (solid line) we show the mean-field result (MF, with the 
pressure determined by quark quasiparticles only) and the mean-field result plus pion and corre- 
sponding scalar contributions. It is evident that at low temperatures the mean-field contribution 
from the quarks is suppressed and the pressure can be described by a free meson gas. Near the 
transition temperature the scalar mesonic modes give a small additional contribution. Finally, 
above temperatures T > 1.5 T c the mesonic contributions become negligible and the quark-gluon 
mean fields dominate the pressure. 

12 The plot shown in Fig. 5 uses the parameters of scenario I. The difference between pressure curves calculated 
from parameter sets I and II turns out to be negligibly small. 
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T/T c 



Figure 6: Trace anomaly (e — 3P)/T A . Dash- 
dotted curve: mean-field results. Solid curve: 
full calculation with inclusion of pions, kaons 
and scalar modes (results for scenario I). 
For comparison, 3-flavor lattice results from 
Ref. [20] are shown. 




T/T c 

Figure 7: Energy density e/T 4 for the mean- 
field and RPA-case. (Legend as in Fig. 6.) 
Lattice data are borrowed from Ref. [20]. 



The RPA treatment of mesonic contributions to the pressure allows us, in addition, to get 
the corrections to the chiral condensates. Exploiting the definition of the chiral condensate, 

m J&A^^e-^ dm q ' ( ' ' 

The pionic corrections to (uu) is computed by differentiating the pion pressure of Eq. (3.11) 
with respect to the up-quark current quark mass: 

dP 

5Auu) = (3.13) 

It turns out, as expected, that the modification of the chiral condensate owing to pions is very 
similar to the results from chiral perturbation theory (cf., e.g., Ref. [45]). At temperatures 
below T c pions tend to soften the condensate and make the chiral transition smoother in the 
range 0.5T C <T <T C . 

Further quantities of interest are the energy density e and the trace anomaly, (e — 3P)/T 4 . 
The trace anomaly, in particular, is relevant here since it is the quantity which can be directly 
computed in lattice simulations (Ref. [20]). This quantity, representing the trace of the energy- 
momentum tensor, is the "interaction" measure which can be expressed in terms of a derivative 
of the pressure with respect to temperature: 

e-3P =T d (P 



T 4 dT \ T 4 

From there it is straightforward to calculate the energy density. Figs. 6 and 7 show the results 
that follow from Fig. 5 both for the mean-field case and with the additional inclusion of mesonic 
(pionic and kaonic) contributions. Our results are compared to 3-flavor lattice data of Ref. [20]. 

3.6 PNJL thermodynamics at finite quark chemical potential 

The nonlocal PNJL approach described in this work can be extended to finite quark chemical 
potential ft = diag(/i M , hj, /j, s ). We do this here with the aim of drawing a schematic phase 
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Figure 8: Chiral order parameter a u for sce- 
nario I, normalized to its value a® at T = ji = 
0, as a function of temperature and u-quark 
chemical potential (note (jl u = /id, fJ- s = 0). 
The thick blue line shows the border for the 
first-order transition. 



Figure 9: Same as in Fig. 8, for param- 
eters of scenario II. The major difference 
between Figs. 8 and 9 is the location 
of the critical point: (a*cep, Tcep)® = 
(180MeV,170MeV) compared to 

(mcep,T CE p) (ii) = (115MeV,195MeV). 



diagram in the (T, fi u ) plane. We set fi = diag(^/, fj,, 0) (i. e. we work in the isospin symmetric 
case n u = [id and with \x s = 0). For the sake of simplicity we restrict ourselves to a scenario 
without diquark condensates. 

The introduction of a chemical potential 13 is accomplished using the prescriptions of the 
Matsubara formalism (see Ref. [43] and Eq. (3.9)): shift the frequencies ujf :n — > ojf n — ifif 
in the particle sector (i. e. in the upper-left submatrix) of the Nambu-Gor'kov propagator 
(3.9) and replace uj^ n — > Uf >n + i^j in the corresponding antiparticle sector (i.e. the lower- 
right submatrix). It is then straightforward to compute the thermodynamic potential £l(T,ft) 
at nonzero ft, following Eqs. (3.8), (3.9) with Matsubara frequencies properly shifted by the 
chemical potential. 

We focus here on the T and [i u dependence of the scalar field a u that acts as a chiral 
order parameter, deduced from the condition -^fe^ = 0. The results are shown in Figs. 8 
and 9 for scenarios I and II, respectively. The profile of a u displays once again the chiral 
crossover transition at [i u = 0. It turns into a first-order phase transition at a critical point 
(located at Tcep ~ 170 MeV and /xcep ~ 180 MeV for scenario I and at Tcep ~ 195 MeV and 
A*cep ~ 110 MeV for scenario II). This qualitative feature is typical for NJL or PNJL type 
models with or without explicit diquark degrees of freedom (see e.g. Refs. [9, 10,48]). Related 
work is reported in Refs. [28, 29] where the nonlocality of the fermionic interaction is introduced 
only in the three-momentum sector. The critical points found in the present calculations differ 
both in their T and [i u values from the aforementioned references. Reasons for such differences 
are given in Refs. [10, 41] where it is pointed out and demonstrated that the location of the 
critical point is extremely sensitive to model details and input parameters. 

The projection of Figs. 8 and 9 onto the T-fi u plane gives the phase diagram of the nonlocal 
PNJL model, Figs. 10. At low \x this phase diagram shows the chiral and deconfinement crossover 
transitions in close contact as already discussed. The deconfinement transition is displayed here 
as a band bounded by dashed lines where the upper and lower bounds are given by values 
$ = 0.5 and $ = 0.3 of the Polyakov loop, respectively, reflecting the relatively soft crossover of 

13 For convenience we use a quark chemical potential throughout this work. The corresponding baryon chemical 
potential is three times the u-quark chemical potential, i. e. fiB = 3/tt u . 
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Figure 10: Phase diagram calculated within the nonlocal PNJL model using the parameters of 
scenario I (left picture) and scenario II (right picture). The solid blue line shows the first-order 
chiral transition (the star denotes the critical end point). The short-dashed blue line marks the 
(chiral) crossover transition while the long-dashed black lines correspond to the deconfinement 
transition (the lower and upper lines correspond to <3? = 0.3 and <!> = 0.5, respectively). 



this transition (see also Fig. 4). At larger values of the chemical potential, beyond the critical 
point, a separation between the chiral and deconfinement transition takes place. 

The area between the (first-order) chiral phase transition and the deconfinement crossover 
has recently been interpreted in terms of a "quarkyonic" phase [50]. It should be pointed out, 
however, that, at nonzero baryon densities and low temperature, PNJL type models are only 
schematic and cannot be considered as realistic. From Fig. 10 it appears that the chiral first- 
order transition boundary meets the /i axis at T = for values of the baryon chemical potential 
as small as \xb = 3/x < 0.9 GeV. This is the domain of nuclear matter that is known to be 
a Fermi liquid of nucleons. The PNJL model works instead with quarks as quasiparticles, the 
"wrong" degrees of freedom in this low-temperature phase at moderate baryon densities. 

The deconfinement transition band at large /i has been calculated using the Polyakov-loop 
effective potential (3.5). This effective potential does not include higher order effects due to 
the presence of quarks at nonzero chemical potential. Such additional /x-dependent effects are 
expected to move the crossover boundary to lower temperatures as \x increases [51, 52]. 

Fig. 10 demonstrates that the position of the critical point is sensitive to small changes of 
the four-fermion coupling G. The relatively small increase of this coupling strength between 
scenarios I and II, keeping hadronic vacuum properties at T = [i = almost unchanged, results 
nevertheless in a significant shift of the critical point in the (T, fi) plane. 

We also confirm that the phase structure is sensitive to variations of the 't Hooft interaction 
coupling strength H (see Fig. 11). An increase of H by only 5% turns the chiral transition 
even at fj, u = into a first-order transition. This behavior might be expected considering the 
Columbia plot (Ref. [53]). On the other hand, a decrease of H, corresponding to a reduced 
rj' mass in the thermal medium, moves the end point to higher chemical potentials and lower 
temperatures. The sensitivity to the axial anomaly observed here in the nonlocal PNJL model 
is, however, less pronounced than that in the local model Ref. [41,54]. We do not observe that 
the end point is removed altogether from the phase diagram as quickly as in the local PNJL 
model. 

Finally, using Eq. (3.8') we calculate the pressure P = — fi at finite chemical potential. In 
Fig. 12 the (normalized) pressure difference 

AP(T, Mu ) := P(T, fx u ) - P(T, fi u = 0) (3.14) 
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is shown for selected values of fi u = pd and compared to (2-flavor) lattice data from Ref. [57]. 
Fig. 13 displays the full result in the T-u u plane. Both figures have been obtained using the 
parameter set of scenario I. 

As already mentioned, the pictures drawn in Figs. 8, 9 and 10 are to be taken as only 
schematic, for several reasons. First, the location of the critical point is sensitive not only to the 
coupling strengths G and H, but also to the input current quark mass [10]. Second, the almost 
constant behavior of a u (T = 0, fi u ) with increasing quark chemical potential is unrealistic in the 
absence of explicit baryon (nucleon) degrees of freedom including their interactions. 

What is actually required as a starting point for extensions to nonzero chemical potential 
is a realistic equation of state at finite baryon density, incorporating the known properties of 
equilibrium and compressed nuclear matter. In such a framework [56], the density dependence 
of the chiral condensate (uu) (or of the scalar field a u ) is well-known to be quite different from 
the profile shown in Figs. 8, 9. The magnitude of (uu) decreases linearly with density p [55], 
with a slope controlled by the pion-nucleon sigma term, and then stabilizes at densities above 
normal nuclear matter through a combination of two- and three-body correlations and Pauli 
blocking effects. The transition towards chiral symmetry restoration at T = is shifted beyond 
at least twice the density of normal nuclear matter [56] . 

Irrespective of these comments, the nonlocal PNJL approach is obviously instructive in 
modeling the chiral and deconfinement thermodynamics at fi = 0. Dealing with finite baryon 
density requires ultimately yet another synthesis, namely a matching of PNJL above the chiral 
transition and in-medium chiral effective field theory with baryons below that transition. 

4 Conclusions and outlook 

We summarize our findings as follows: 

i) The nonlocal generalization of the PNJL model to iVf = 3 flavors incorporates all important 
nonperturbative features of low-energy QCD with inclusion of strange quarks: spontaneous 
and explicit chiral symmetry breaking, the axial anomaly, and thermodynamical aspects of 
confinement in terms of the Polyakov loop. A separable form of the underlying nonlocal effective 
interactions between quarks, including the axial U(1)a breaking six-fermion vertex, proves to 
be successful in reproducing the physics of the pseudoscalar meson nonet. Chiral low-energy 
theorems are shown to be fulfilled. 

ii) At the level of generalized gap equations, the resulting momentum-dependent quark 
(quasiparticle) masses permit establishing connections with instanton physics, lattice QCD and 
(Landau gauge) Schwinger-Dyson approaches. The nonlocality distribution of the effective inter- 
action between quarks reflects typical instanton sizes of about 1/3 fm. This distribution replaces 
the artificial sharp momentum space cutoff in standard (local) NJL type models and links the 
quark mass function M(p) smoothly to the correct QCD behavior at large momentum scales. 

iii) The thermodynamics of the three-flavor nonlocal PNJL model reproduces corresponding 
Nf = 2 + 1 lattice QCD results with almost physical quark masses surprisingly well. In particu- 
lar, the dynamical entanglement of the chiral and deconfinement crossover transitions observed 
previously in two-flavor PNJL models is confirmed also for Nf = 3. The interaction measure 
e — 3P is well reproduced around the critical temperature T c ~ 0.2 GeV. Mesonic contributions 
to the pressure can be systematically incorporated. 

iv) The critical point in the phase diagram and its location in the plane of temperature T 
and baryon chemical potential hb remain an open issue. The position of this critical point turns 
out to be extremely sensitive to fine-tunings of parameters (quark masses, coupling constants), 
as already found in previous studies. 

v) Previous calculations at large quark chemical potentials, using local PNJL models, were 
limited by the momentum space cutoffs characteristic of such models. While these limitations 
are overcome in the nonlocal PNJL approach, one must still be aware of the fact that, at low 
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Figure 11: Location of the critical point de- 
pending on several values of the 't Hooft cou- 
pling strength H in units of the coupling 
strength Hq of scenario II. 



Figure 12: Pressure difference at finite 
chemical u-quark potential, AP(T,fi u ) = 
P(T, u, u ) — P(T, 0) (parameters of scenario I) 
for selected values of \i u compared to lattice 
data from Ref. [57]. 




Figure 13: Excess pressure AP(T,fi u ) = P(T,fj, u ) — P(T, 0) at finite chemical u-quark potential 
calculated in the nonlocal PNJL model in units of T 4 for scenario I. The thick solid blue line 
indicates the chiral crossover and first-order transition line. 



temperatures and moderate baryon densities, PNJL models do not operate with the proper 
nucleon degrees of freedom relevant at such densities. In the baryonic phase around and below 
~ 1 GeV, the phase diagram should represent a nuclear Fermi liquid and not a quarkyonic 
quasiparticle system. 

Implementing the constraints from realistic nuclear or neutron matter equations of state on 
the QCD phase diagram remains a challenge for the future. 
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A Derivation of the 't Hooft interaction 



In this appendix we show how the Kobayashi-Maskawa-'t Hooft determinant expression, Eq. (2.6), 
can be cast into the form used in this work, Eq. (2.8). 

In order to write the 't Hooft determinant in a more tractable way, we apply Newton's and 
Girard's formula 



det J* = I (tr^) 3 - \ (trj±) (W ±2 ) + 



(A.l) 



Here tr indicates the trace over flavor space only. We use the Gell-Mann matrices as a basis in 
flavor space, {Ao, Ai, . . . , Xs}, with the additional definition Ao := \J\ diag(l, 1, 1) in order to 
maintain tr{A Q • Xp} = 25 a f3 for all a, (3 e {0, ... ,8}. This allows us to write 



J^ 1 — C t^a 



tllXaJ^ = 2C- 



-a 



(A.2) 



a=0 



and, consequently, = ^^{XaJ^}. 



Furthermore, from Eq. (2.6) we have with the definitions (2.5) 



-trlA^z)} { 



' / d 4 zA^(x + |)(lT75)C(z)^( 



X 



By means of Eq. (A.2), this allows to write = ± ^j p , or, inversely jf = 2 (c+ + c~) ,j p = 
-2i(c+-c-). 

Next, we return to Newton's and Girard's formula, Eq. (A.l), and use 

tr (j ±2 ) = tr (ctX a c±Xp) = <£<%26 af) = 2c±c± 
tr (j ±3 ^j = c t c^c^tr (A Q ,A /3 A 7 ) . 
Inserting this into Eq. (A.l), one has 



det J + + det J~ 



-4 



+3 _3 
Cq -+■ Cq 



+ 



C a Cg c 7 + c a c 7 



tr (A a A,gA 7 ) . 



Now, using the relations between the c's and the currents j s ,j p , we have 



„+ 3 



2Re 



.+ 3 



c C a c a + c C a C o — ^Re [c^C+C+j — ^ 



^io 1 (io ~~ 3jo 
l 



JO [J a J a J Z J0 JaJa 



C a Cg c 7 + C a Cg c 7 



2Re 



c a c j3 C 7 



Inserting this in the previous formula leads to 



1 



;P 2 



16 V 2 



doU + + riot J" --\/^fif-3ir)-^^ 



+ 
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where summation over a, (3, 7 G {0, . . . , 8} is implicit. 

Finally, using the SU(3) structure constants f k e m , d k e m , defined through [X k , A^] = 2if k g m X m 
and {X k , X e } = ^S ki + 2d Mrn X m , respectively, one obtains X k X e = i f k £ m X m + d kim X m + \5 ki and, 
hence, 

tr (A fc A^Ai) = 2 i f k £ m 5 mi + 2d k i m 5 mi (A. 3) 

(for k, £, m, % G {1, . . . , 8}) which allows us to write 

det J + + det J- = -L ^1 (jf - 3i 5 i p2 ) - ^ ^| ( j 5 jf £ - j!££ - 2j s k jo£) + 
+ -^d ekrn (jfj k jm - 3jfj k jm) , 

where, again, £, k, m G {1, . . . , 8}. 
If one sets 

A*/3 7 := y £ijk£mnt {K) irn {>>p) jn (\) kl for a, (3, 7 G {0, . . . , 8}, (A.4) 
then the expression above can be written in a more compact form as 

det J + + det J" = j^A a ^ - 3£j£jP) . (A.5) 

This is the form given in Eq. (2.8), with the coupling constant properly adjusted. 

B Second-order contributions to the action 

We demonstrate how the expressions for the inverse meson propagators, Gp, can be derived. 
The self-energy contribution IJij can easily be derived from the fermion determinant, Eq. (2.12'), 
using the following formulas for the functional derivatives 

6 dn det J = -A-Tr In J = Tr J^" W^'^ ^ 



5iTij(k) STTij(k) \ ' SiTij(k) 

5 , , ,6 J 



dirij(k) Stt, 



Moreover, owing to the SPA equations, Eq. (2.13), the auxiliary fields 14 S a ,P a are implicit 
functions of a and ir. This implies calculating the second derivative of the expression 



G H 



•Se '■— <?aS a + TTa-Pa + ~^ (S a S a + P a P a ) + — A a /3-y (S a SfsS y — 3S a PfiP, 



1) 



Neglecting first the space dependence of the fields we may first introduce the matrices 
■^a a X a and V' = ^vr a A a . The SPA equations (2.13) in this new basis then read 



3H 

+ GS a X a + — AafrXaiSpS^ - PpP 7 ) = (B.la) 
yW + GP a X a - ^-A a ^X a S^ = 0. (B.lb) 



From the first derivative of Eq. (B.la), 



4 from now on we omit the tildes on S a , P a 
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it follows that = for all a G {0, ... , 8} and i,j G {1, 2, 3}, recalling that P a = for all a 

in mean-field approximation. 

The second derivative of Eq. (B.la) leads to 

r} PSa m 5*S, _ 3H 5P P 5P, 

LrAa l r -Tr^ap-yAaap- — -— A. a ^A a - . K^- 1 ) 

0TT k e OTTij 2 0TT k l dlTij 2 dlT k g OTTij 

Analogously, one has from the second equation 

V26 im 6 jn + {K) mn - ^-A aPl {K) mn S r 6 ' - " 



or, by contraction with (A e ) r 



5P e 3H 5P, _ 1 

" — A ^ s ^ ~ -71 (Ae) « • (R3) 

Finally, from the second derivative 

and it follows that S J' = for all a, i,j, k,£ in mean-field approximation. 
The sum of the SPA equations gives 

3H 

cr a S a + n a P a + G(S a S a + P a P a ) + -^A a p 1 (S a SpS 1 — 3S a PpPj) = , 
so that one can write 

1 H 

Se = —-G(S a S a + P a P a ) ^-Aap^^aSpS-f — SSaPpPy) . 

Finally, applying identities (B.2) and (B.3) we may deduce the desired derivative 

5 2 S E 1 n , SPp 



We conclude that the additional term is given by the solution of relation (B.3) contracted by 

Wmn 2~A*/3 7 ^ a ">mn S P~fa~ L = -^^im^jn 5 

this can be further simplified by noting S a = ^tr(\ a S) and 

6P 1 1 5 Pry 

A*/3 7 ( A a)m„ ^PT~ = ^ e rsk£uv£ {K) ru (^) sv (A 7 ) M ~S t (\p) tt (A Q ) 



1 ^ ( \ \ C 

^ntk^mti-^- \^y)k£ <->t ■ 



Consequently, the equation to be solved is 



Defining (r ij>mn ) 1 := ^ {K) mn we may write 

8 2 Se 



29 



where solves the system given in Eq. (2.28). 

Finally, we consider the functional derivative of terms of the form f d x S a (x)Sp(x)Sj(x) 
etc. The first derivative with respect to ifij{y) generates a 5 function, S(x — y), hence 

J d 4 xS a (x)^(x)S 7 (x) -> S a (y)Sf,(y)S y (y) . 

The second derivative with respect to vr^(z) generates an additional 5(y — z). This means that in 
mean-field approximation the functional dependence of the fields after a Fourier transformation 
is given by 

Treating analogously the contributions from the a field, we arrive at Eq. (2.25). 
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